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We write down and solve a closed set of Schwinger-Dyson equations for the two- 
matrix model in the large iV limit. Our elementary method yields exact solutions for 
correlation functions involving angular degrees of freedom whose calculation was impossible 
with previously known techniques. The result sustains the hope that more complicated 
matrix models important for lattice string theory and QCD may also be solvable despite the 
problem of the angular integrations. As an application of our method we briefly discuss 
the calculation of wavefunctions with general matter boundary conditions for the Ising 
model coupled to 2D quantum gravity. Some novel insights into the relationship between 
lattice and continuum boundary conditions are obtained. 
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Matrix models in the large N limit are widely believed to be capable of giving qual- 
itative and even quantitative insight into non-abelian gauge theories. The simplest such 
models provided early on very explicit solutions to low-dimensional QCD Over the 

past few years another physical application of large N matrix models has emerged: It 
has become clear following that hermitian models give a precise covariantly regularized 
definition of perturbative bosonic string theory and 2D quantum gravity. For the first 
time one was able to attempt to address even non-perturbative issues in string theory Q . 
However the class of solved cases remains too small: The dimension of the target spaces in 
which the world sheet (made up from the matrix model Feynman diagrams) is embedded 
is too low. In both the QCD as well as the string theory case there is simply not enough 
room for the physically most interesting excitations. It is well known that the the techni- 
cal obstacle preventing further progress consists in our inability to deal with the relative 
"angles" between interacting matrices in physically interesting multi-matrix models. This 
inability to rid oneself of the angles would already be true in the simplest interacting case 
— the two-matrix model — were it not for the Itzykson-Zuber formula enabling one to 
actually integrate them out || . Unfortunately it appears impossible to generalize this ap- 
proach to any target space lattice containing even one cycle. The question remains whether 
there does not exist a general approach which does not break down once the angles can no 
longer be eliminated. In principle such a method is known: The loop (Schwinger-Dyson) 
equations ||. It is therefore of some interest to investigate whether the two-matrix model 
can be understood in this language. It will be demonstrated below that this is indeed 
the case. After deriving a set of closed loop equations using elementary combinatorial 
reasoning we derive the algebraic (for polynomial potentials) equation for the resolvent. 
This resolvent is called "masterfield" below since it encodes, via a Hilbert-transform, the 
information about the spectral density of the matrices. We then show how the general one 
trace correlators are expressed as rational functions of the master field (further justifying 
the use of language). We finally apply our formulae to 2D quantum gravity and sketch 
how to extract some hitherto unavailable information. 
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Consider then the two-matrix integral 



Z = jvAVB e -A f T r [V(A)+V(B)-cAB] (1) 

with A, B hermitian in the large iV limit. V(A) is a polynomial potential which we will 
take for the most partB to be V(A) = ^A 2 — ^A 3 . The simplest "observables" in this 
model are the moments 

W n = i ( Tr A n ) (2) 

where the expectation values (...) are taken with respect to the ensemble ([I]). These mo- 
ments are the correlators calculable with the standard technique of orthogonal polynomials 
0. This is not the case for more complicated operator expectation values 

Wg% Xi ... tnu , mh = ^ ( Tr A ni B mi . . . A^B- ) (3) 

After changing variables to the eigenvalues of A and B a highly nontrivial angular inte- 
gration remains to be performed for each correlator. Here we will find ©,@ through a 
Schwinger-Dyson technique!!. It is useful to introduce the resolvent 

1 1 oo 

W{P) = ^(Tr ¥ - J ) = J2 (p) n+1 W n (4) 

n=0 



as well as the generating functions 

w^(a, Qi , •-i*.q») = £e v 1 _ Aqi 1 _ B - pTaqTb > 

00 1 1 11 

_ ( , __ N | n i + 1 ('_!_'| m i + 1 ( _)ri k + l ( ^_\m k + l w(2fc) 

rii,mi=0 " 



(5) 



1 Everything works equally well for more general polynomials and for the case where A and 
B are controlled by non-equal potentials. Some remarks on the first generalization follow below. 

2 Our method could be extended to the most general connected correlators involving an arbi- 
trary number of such traces as well as to include ™ corrections. This will not be done at present 
since our main point is merely to demonstrate how to circumvent the "angle" problem. 



2 



and the auxiliary functions 

1 1 oo 

W^P) = - ( Tr ¥ - x ) = J2 (p) n+1 W?J (6) 

n=0 

The strategy is now to write down recursion relations between the coefficients (|2]), ([|) and 
use them to find equations relating the various generating functions. The recursion for W n 
is (n > 1) 

n-2 

W n =g W n+1 + c W™ ltl + Wj (7) 

J=0 

Here c, g are defined through eq.(|l|) above. Such recursions may be easily visualized and 
derived, using elementary combinatorial arguments, by representing the W' 2k ' by "wheels" 
with 2k alternating sequences of n\ spokes of color A, n<i spokes of color B, and so on. 
Cutting out one spoke and considering the possible new configurations gives the wanted 
relations. The method is explained e.g. in ||. Eq.([7|) implies, using the equation 

(P-g P 2 ) W(P) = 1 - g {W x + P) + W(P) 2 + c W^P) (8) 

For c = (i.e. the one-matrix model) this is the well known basic loop equation 0. In 
the two-matrix case the equation does not close and we need to consider also two further 
recursions: 

<i = 9 W<$ + c W n+1 (9) 

and 

n-2 

W% = 9 W$ ltl + c wi% + ]T W 3 W™ 2 _ jtX (10) 

3=0 

which are valid for n > and n > 1, respectively. Eqs. (|9|) , (|10|) are summed up to give 

W 1 (P)=gW 2 (P)+cPW(P)-c (11) 

and 

(P - gP 2 ) W^P) = (1 - gP) W 1 - gW[ 2 l + cW 2 (P) + W(P) W^P) (12) 
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Now eqs.(|^),(|T^),(|T2|) constitute a closed set of expressions for the unknown functions 
W(P), Wi(P), W2(P) and we may deduce an algebraic equation for the master field: 

W(Pf + \- - 2(P - gP 2 )} W(P) 2 + 
9 

3 

+ [l - g(W x +P) + —P--(P- gP 2 ) + (P - gP 2 ) 2 ] W(P)+ ( 13 ) 

9 9 

+ [(gW 1 +gP-l)(P- gP 2 ) + r x P + r ] = 



It was found by different methods in || , |T0| . As in the case of the one matrix model master 
field equation the remaining constants! Wi,ri,rg are fixed by requiring that W{P) possess 
only one cut on the physical sheet. 

Now let us extend these simple combinatorial arguments to derive some novel results. 
The recursions (|9|) , fllCf) are easily generalized (for n > 1 and all m > 0) to 



n-2 



W£l = g W {2 ^ m + c Wi\ m+1 + J2 Wi Wi%_ j>m (14) 

3=0 

Using eq.(||) one is led after some algebra to the equation 

W {,Q) ~ P-gPi-cQ-W(P) (15j 

Thus, in view of eq.(f§), the first nontrivial correlator involving "angular" degrees of freedom 
is seen to be a simple rational function of the master field W(P). 

Consider the next correlator W^ 4 -*. The basic procedure is identical to the one above, 
but the combinatorics gets a bit more involved. The recursion is (ni > 1) 

w w =q w {4 L + c w {A) , ^ + 

ni ,mi ,ri2 ,m,2 i> m +l,mj ,ri2 ,7112 1 ni — I,mi + l,n2,rri2 1 

ni— 2 n 2 -l (IfW 

+ V w w {4) + V w {2) w {2) 

1 / j "J ni -2- j, mi, n 2 ,»n 2 1 / j j,mi m 2 ,7ii +n 2 — 2— j 
j=0 j=0 



One has n = -c + gWi - g W 2 and r = -g + - - (1 + c)Wi + 2gW 2 - g ' W 3 and we 



have expressed W{ { through W\, r\,ro. 



In order to be able to sum up this recursion and find an equation for one needs a 

further relation (valid for all mi, n 2 , m 2 > 0) 



C W (4) =W {2) -aW {2) 

mi — 1 T712 — 1 

(2) 



- v w j w {2) . , . - y w j w (2) , , . 

/ j 3 n 2 ,m 1 +m,2 — 1— J / j 3 n2,m 1 +m 2 — l—3 



(17) 



3=0 



Then, using properties like Wq 
amount of algebra 

W^(P 1 ,Q 1 ,P 2 ,Q 2 ) 



(4) 



r(2) 



,m 1 ,n 2 ,m 2 W^^+m,, one derives after a considerable 



P x - gP* - cQ 1 - W(P X ) 
Dq(QuQ 2 )- [^{Q-gQ 2 -W(Q 1 )-W(Q 2 ))~l + gP 1 } W^(P 2 ,Q)+ (18) 

+ (c- W^(P 2 ,Q 1 )) D P (P U P 2 ) ■ W^(Q 2 ,P) + 9 - W(P 2 ) 

where we have made use of the "combinatorial derivatives" 

f(P 2 ) - /(Pi) 



D P (P 1 ,P 2 )-f(P) = 



(19) 



P2-P1 

It is now straightforward if somewhat tedious to generalize the recursions (|l6|), (|l7|) and 
thus compute all remaining (k > 3) one trace correlators W^ 2k ^: 

1 



W^(P 1 ,Q 1 ,...,P kl Q h 



■ x 



P 1 -gP 2 -cQ 1 -W(P 1 ) 
D Q (Q u Q k ) ■ {- c {Q-gQ 2 - W{Qx) - W(Q k )) - 1 + gPx]- 

■W^ 2k - 2 \P 2 ,Q 2 ,...,P k ,Q)- 

- 9 -D P {P 2l P k ) ■ W^ 2k ~ 4 \P, Q 2 , . . . , Q fc _x) + 
c 

k-i 

+ 9 - J2(D Q (Q l ,Q l ).W^ 2l - 2 \P 2 , . . . , P h Q))- 

1=2 

■ (D Q (Qi,Q k ) ■ W^ 2k - 2 ^ (P l+1 , ...,P k , Q)) + 

+ cD P (P l ,P 2 ) ■ W^ 2k ~ 2 \P, Q 2 , ■ ■ ■ , Qk)- 

k 

- ^ (2Z " 2) (Pi,Qi,---, Qi-i)D P (P 1 ,P l ) ■ W ^ 2k+2 - 2l \P, Q h ...,Q k ) 



(20) 



We have thus succeeded in recursively expressing all higher one trace correlators as rational 
functions of the master field W(P) alone. Note that the obtained formulae hide the fact 
that the possess a cyclic symmetry with respect to the variables Pi,Q{. 

We now briefly sketch which modifications occur for a general polynomial potential 
of order m + 1: V(A) = — Y1T=2 ~9rA r (here #2 = — !)• This generalization is of interest 
since any minimal C < 1 conformal theory may be obtained by an appropriate!! choice of 
V(A) ||11|| , ||12|| . Studying the necessary recursions (their number increasing with the order 
of the polynomial) one finds the masterfield equation 

W(P) m+1 + hm(P) W(P) m + ... + hx(P) W(P) + ho(P) = (21) 

The hi(P) are known polynomials in P, e.g. h m (P) = — ( c g 9m 1 +mV'(P)). It is interesting 
to observe that the loop functions W(P) are algebraic, just as in the case of the 0(n) model 
at rational points |TjJ. One again finds that the higher one-trace correlators are rational 



functions of the master field; e.g. reads 

W (2 )(P c) _ - Eg* £«=i fr+i P r ~ s WsMQ) - C W(P) 

W tPQ)- V'{P)-cQ-W{P) [ZZ) 

The functions Wi(P) are related to each other and W(P) through expressions similar to 
(Bp; these relations are derived using the same methods as above. Alternatively one can 
find them by investigating ( |2"2"| ) for Q — > oo and using (g). 

Let us complement this analysis with a brief application of our results to continuum 
2D quantum gravity. It is well known that the two-matrix theory eq.(P describes at 
its critical point the C = \ conformal theory coupled to 2D gravity []l4j|,||15|l. Indeed, 
analyzing the master field equation ([13]) one reproduces the correct critical valuesS for 
c,g,P Hl5fl . Scaling g = g* — a fx and P = P* + a z where a is the cutoff and \i and z are the 



4 For non-unitary minimal models one actually needs unequal potentials for the two matrices. 
We refrain from giving formulae for that case. 

5 They are c* = 0.1589 g* = 0.2004 . . ., P* = 3.6850 .... 



continuum bulk and boundary cosmological constants one reproduces the correct tree-level 
one-boundary amplitude |p.6| ,|]l7l, ||18|l (in the first reference this result was obtained from 



a different formulation of the C = \ theory and in the latter references from the same 
model, but in Laplace-transformed form): 

W(P) ~ J w(z) = J [(z + \/z 2 -n)% + (z- Vz 2 - n)*] (23) 

Our calculation of the general moments (|3|) enables one to now study a much richer set of 
matter boundary conditions. However, even though eqs . (P~5|) , QI£|) , (|2T)|) are rather explicit 
rational functions of the amplitude W(P) the scaling of these expressions is more subtle 
than one might expect at first. In particular, the continuum limits of the amplitudes 
are not simply rational functions of the scaled loop function w(z). Instead one has 
to expand W(P) to higher powers in the cutoff a than order |; i.e. one needs to keep terms 
that were neglected in (|23|). Here we will only present results for fi = 0; one finds, scaling! 
Pi = P* + az, Qi = P* + az 

W 2k (P u Q l7 ...,P k ,Q k ) ~ (az) 7 / 3 - k (24) 

It is instructive to compare this scaling law to the continuum KPZ/DDK prediction for 
the insertion of 2k boundary operators Oi into a disc ((p is the Liouville field, 7 = \/l and 



a t the dressing charge of O t ) Pf,p0|,pi 

2k 



f[ (fO,e a ^)\ ~ ^-fE-l^ (25) 



i=i 

Therefore one may deduce the charge of the continuum boundary operator O producing 
the behavior @D to be ^ = \. One then finds, using ^ = §(7 -VI + 48 A) [0,0, [|T[, 
the flat space boundary scaling dimension A of O: 

A(O) = ^ (26) 



6 Note that one could also weigh each boundary segment with its own boundary cosmological 
constant. 



This is not a dimension occuring in the minimal Kac table of the Ising model. It ap- 
pears however in the extended Kac table, suggesting that it does correspond to a nonlocal 



suit might be surprising at first: An operator (TrA n B m ) would naively correspond to a 
boundary with a sequence of up spins followed by a sequence of down spins. It is well 
known that such a spin configuration corresponds to the insertion of two spin operators 



(p6|). We believe that the key to resolving this puzzle lies in considering more carefully 
the original random surface interpretation of the theory (|l|). For closed surfaces we can 
imagine the Ising spins to sit either on the vertices of the </? 3 -lattice or alternatively on 
the vertices of the triangulation dual to the </? 3 -lattice: Kramers- Wannier duality holds 
also for random lattices. Once we include boundaries it is easily seen that insertion of the 
operator Tr(A n + B n ) into ([!]) corresponds to free boundary conditions for the spins on 
the triangulation. An operator like (TrA n B m ), on the other hand, does not correspond 
to a simple local boundary condition for the Ising spins on the triangulation. In short, 
we would like to argue that Kramers- Wannier duality is destroyed once we consider open 
surfaces and that we should imagine the spins to sit on the lattice dual to the ip 3 lattice. 
Further evidence for this interpretation may be obtained from the loop gas approach || 
which enables one to consider boundaries with 2k alternating segments of up and down 
spins (situated at the vertices of the triangular lattice) separated by domain walls. A sim- 
ple analysis |23| of the disc endowed with these boundary conditions along the lines of ||24|| 
gives the behavior ~ zi~i k . This agrees (using A(cr) = h) with the continuum prediction 



Obviously further results (e.g. for \i ^ and for the higher critical points) could be derived 
from the above formulae. 

We have thus demonstrated that the master field of the two-matrix model with poly- 
nomial potential is an algebraic function and that the one-trace correlators are rational 



boundary operator in a non-minimal extension of the C = \ conformal theory. The re 



a into the boundary j52|. But the boundary dimension of a is A(cr) 



\ as opposed to 




(27) 
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functions of the master field. It is not easy to imagine how to obtain the latter result 
using the standard method of orthogonal polynomials: Some rather complicated angular 
integrations have to be performed (but see ||25|| ). The most interesting question is of course 
whether the study of the loop equations of more complicated matrix models will lead to 
new solutions. It would be of some interest to understand for which target space lattices 
the master field, be it for hermitian or unitary matrices^, remains algebraic. There are 
rather "small" matrix models like hermitean five matrix models and unitary Eguchi-Kawai 
type models |§ whose solution would presumably lead to important new physical insights: 
The former would contain a C > 1 critical point in bosonic string theory and the latter 
would give an analytic expression for strongly coupled QCD. Is it truly impossible to find 
the master field for these models? 
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7 We remind the reader that for sufficiently simple link structures unitary models possess an 
algebraic spectral distribution as well @,[26]. Of course we know that already a model consisting 
of four linearly interacting hermitean matrices can not possibly have an algebraic master field: 
Corresponding to a C = 1 lattice theory (the 4-state Potts model), it will certainly have a 
logarithmic singularity in its eigenvalue density. 

9 



References 



[1] E. Brezin, C. Itzykson, G. Parisi, and J.-B. Zuber, Commun. math. Phys. 59 (1978) 
35. 

[2] D. Gross and E. Witten, Phys. Rev. D21 (1980) 446. 

[3] V. Kazakov, Phys. Lett. 150B (1985) 282; F. David, Nucl. Phys. B 257 (1985) 45; V. 

Kazakov, I. Rostov and A.A. Migdal, Phys. Lett. 157B (1985) 295; J. Ambjorn, B. 

Durhuus, and J. Frohlich, Nucl. Phys. B 257 (1985) 433. 
[4] E. Brezin and V. Razakov, Phys. Lett. 236B (1990) 144; M. Douglas and S. Shenker, 

Nucl. Phys. B335 (1990) 635; D. Gross and A. Migdal, Phys. Rev. Lett. 64 (1990) 127. 
[5] C. Itzykson and J.-B. Zuber, J. Math. Phys. 21 (1980) 411. 
[6] For a review see: A.A. Migdal, Phys. Reports 102 (4) (1983). 
[7] M.L. Mehta, Comm. Math. Phys. 79 (1981) 327. 
[8] I.R. Rostov, Nucl. Phys. B326 (1989) 583. 

[9] J. Alfaro, "The large N limit of the two hermitean matrix model by the hidden BRST 
method", CERN preprint CERN-TH-6531-92, |hep-th/9207097| . 
[10] D.V. Boulatov, "Infinite tension strings at D > 1", N. Bohr Institute preprint NBI- 
HE-92-78, |hep-th/ 92 11064 . 



[11] M.Douglas, "The two-matrix model", 1990 Cargese workshop; E.Martinec, Comm. 

Math. Phys. 138 (1991) 437; T.Tada, Phys.Lett. B259 (1991) 442. 
[12] J.-M. Daul, V.A. Razakov and I.R. Rostov, in preparation. 
[13] I.R. Rostov and M. Staudacher, Nucl. Phys. B384 (1992) 459. 
[14] V.A. Razakov, Phys. Lett. 119A (1986) 140. 
[15] D.V. Boulatov and V.A. Razakov, Phys. Lett. 186B (1987) 379. 
[16] I.R. Rostov, Phys. Lett. 266B (1991) 317. 
[17] M. Goulian, Phys. Lett. 264B (1991) 292. 

[18] G. Moore, N. Seiberg and M. Staudacher, Nucl. Phys. B 362 (1991) 665. 
[19] V. Rnizhnik, A. Polyakov and A. Zamolodchikov, Mod. Phys. Lett. A3 (1988) 819. 
[20] F. David, Mod. Phys. Lett. A3 (1988) 1651; J. Distler and H. Rawai, Nucl. Phys. 
B321 (1989) 509. 

[21] E. Martinec, G. Moore and N. Seiberg, Phys. Lett. B263 (1991) 190. 
[22] J. Cardy, Nucl.Phys. B324 (1989) 581. 
[23] M. Staudacher, unpublished. 

[24] B. Duplantier and I. Rostov, Nucl. Phys. B340 (1990) 491. 

[25] S.L. Shatashvili, "Correlation functions in the Itzykson-Zuber model", Princeton-IAS 

preprint IASSNS-HEP-92-61, [hep-th/9209083| . 
[26] D. Friedan, Commun. math. Phys. 78 (1981) 353. 



10 



